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, We discuss mechanisms of decoherence of a confined Bose-Einstein condensate at finite temper- 

■ atures under the explicit condition of conservation of the total number of bosons in the trap. A 

criterion for the irreversible decay of the condensate two-time correlator is formulated in terms of 
the Orthogonality Catastrophe (OC) for the exact N-body eigenstates, so that no irreversible decay 
^ l' occurs without the OC. If an infinite external bath contacts a finite condensate, the OC should 

D ' practically always occur as long as the bath degrees of freedom are interacting with each other. We 

find that, if no external bath is present and the role of the bath is played by the normal component, 
\^ ' no irreversible decay occurs. We discuss the role of the effect of level repulsion in eliminating the 

OC. At finite temperatures, the time-correlations of the condensate isolated from the environment 
are dominated by reversible dephasing which results from the thermal ensemble averaging over real- 
^^-H ' izations of the normal component. Accordingly, the correlator exhibits gaussian decay with certain 

O , decay time Td which depends on temperature as well as on intensity of the shot noise determined 

by the statistical uncertainty in the number of bosons N in the trap. 
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^ . I. INTRODUCTION 

o 



Creation of trapped atomic Bose-Einstein condensates has made possible experimental study of the fundamental 
concepts which so far were tested mostly in "thought experiments". Many intriguing questions are associated with 
the finiteness of the number of bosons N forming the condensate. As discussed by Leggett and Sols , depending 
[ on the environment, the Josephson phase of a finite capacity Josephson junction may exhibit either a ballistic or a 
^ I ■ diffusive spread in time. For example, if the condensate randomly exchanges bosons with some environment, phase 

■ diffusion becomes completely irreversible 
— . , A dephasing of the wave function of the confined condensate - the so called phase diffusion - has been discussed 
in Refs. |^-|^ in recent times. The nature of this effect is closely related to the concept of broken gauge symmetry 
, associated with the formation of the condensate wave function {'^{x,t)) ^ 0, where 'i/{x,t) is the Bose field operator, 
' and the averaging is performed over the initial state taken as the coherent state. At zero temperature, the phase 
' diffusion refers to the decay of {9{x,t)) due to the two-body interaction [Q-^. This effect is purely quantum, and 
^ ■ occurs without any bath. Such a decay was predicted to be reversible in time, so that spontaneous revivals of the 
S ' wave function should occur The phase diffusion effect in an atomic condensate which is completely separated 

. from the environment was considered by Graham |0J|], and it has been suggested that {'^{x,t)) decays irreversibly 
(-H ■ due to the particle exchange between the condensate and the normal component confined in the same trap. 
Q ', The experimental study of the temporal phase correlations has been conducted by the JILA group . It has been 
O found that no detectable decay of the phase exists on the time scale of the experiment < 100ms. The question was 
^ _ then posed in |^ why the phase correlations are so robust despite an apparent fast relaxation of other degrees of 
freedom such as, e.g., the relative motion of the two condensates. 

In the present work, we address the question P] of the phase correlations. However, this question is not analyzed 
^ , in the context of the problem of decoherence of [^{x,t)). Strictly speaking, as long as the operator 'i'{x,t) describes 
conserved particles, the mean {9{x,t)) is not physically observable. We analyze a physically measurable quantity 
which is the two-time (and space) correlator 

p(x,x',i,t') = (vl/t(x,i)*(x',0). (1.1) 
This correlator can be m easu red by, e.g., scattering of fast atoms PQ]. We formulate a condition for the irreversible 



decay of the correlator (1.1) in terms of the Orthogonality Catastrophe (OC) for the projected N-body eigenstates 
(defined below). 

We consider a situation when the one-particle density matrix (OPDM) p(x, x', i) = p(x, x', i, t) exhibits equilibrium 
off-diagonal long range order (ODLRO) Q ( see also the discussion in Ref. on how the ODLRO should be def ined 



in the atomic traps). We ask how the existence of the ODLRO affects the temporal behavior of the correlator (1.1) 
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as t — t' ^ +00. It is worth stressing that this question should not be identified with the problem of the decay of 
the broken U{1) symmetry state, because such a state is not a prerequisite for ODLRO |jl^. This difference becomes 
especially important while studying the long time correlations in large (but finite) N condensates. As was found 
by Wright et al. Q|, the interaction between bosons introduces some break time ^ ^/N ||T^, so that employing the 
classical field (\E'(x, i)) as a proper variable obeying the Gross-Pitaevskii equation at times longer than this break 
time becomes no longer valid. A similar situation occurs in the description of the normal excitations, when the mean 
field approach breaks down at long times , with the break time being ^ y/N as well. The limitations of the 

Gross-Pitaevskii equation at long times have also been pointed out by Castin and Dum ||l^ . If one is interested in 
the long time evolution, the correlator (1.1) must be analyzed. 

In this paper we are interested in the long time correlati ons which, however, can be quite measurable in the current 
traps. We calculate the condensate part of the correlator (^]^) in equ ilibrium at T 7^ and show that, as long as no 
external bath is present, the time dependence of the correlator ( |l . l[) can be represented as a result of the thermal 
averaging of the non — decaying exponents where each one can be viewed as representing a specific realization of the 
chemical potential (defined below) for the exact A'^-body eigenenergies. This averaging results in the dephasing of 
the correlator. The time Td of this dephasing is determined by the thermal fluctuations of the normal component for 
fixed N. In any realistic experiment, the value of N varies statistically from run to run, and produces shot noise, also 
contributing to the dephasing rate ■ The equilibrium condensate evolution can be viewed as though it occurs in 
the frozen environment created by the normal component. The dephasing, then, results from the ensemble averaging 
over the possible realizations of the normal component and N. 

In the following we will formulate a necessary condition for the dissipative decay of the correlator (1.1) in terms of 
the Orthogonality Catastrophe for the exact A^-body eigenstates. 



II. ORTHOGONALITY CATASTROPHE AS A NECESSARY CONDITION FOR THE DISSIPATIVE 
DECAY OF THE CONDENSATE TWO-TIME CORRELATOR AT T 7^ 



The correlator (1.1) carries the most complete information about confined bosons. The averaging (...) in Eq.( |l.l| ) 



is performed over the initial state (or states). In what follows we will assume that the system of N bosons is in 
the thermal state, so that the averaging is performed over the thermal ensemble of the exact iV-body eigenstates 
|m, N) and the eigenenergies Em{N), where m denotes a set of the quantum numbers specifying the eigenstate with 
given N of the many body Hamiltonian H. In other words, if the ground state |0,A^) corresponds to the pure 
condensate (characterized by the energy Eo{N)), the states with m 7^ describe the condensate and the normal 
component in the trap, with m carrying the meaning of the set of all possible quantum numbers of the normal 
excitations. Should some external bath be present in addition to the normal component, providing (coining the 
terminology Q) the intr insi c bath, these numbers characterize both the normal component and the bath. In the 
thermal equilibrium Eq. (|l.l[ ) is reduced to p(x, x',i, i') = p{x.,x',t — t',0). Accordingly, the OFDM becomes time 
independent p(x, x',i) = /9(x, x', 0, 0). In order to simplify notation, we will set t' = 0. 

As a matter of fact, for temperatures T below and not very close to the Bose-condensation temperature Tc, the 
behavior of ( |l.lD is dominated by the condensate part. The condensate part po(x, x', t, 0) of p(x, x', t, 0) is defined in 
terms of some macroscopically populated eigenstate (po (x) of the OFDM ||ll|,|l^ . In the case of the weakly interacting 
Bose gas, /9o(x, x', 0) can be selected by employing the following standard representation 

*(x,t) =ao(0¥'o(x) +?A'(x,0, (2.1) 

where ao{t) removes one boson from the condensate, and ip' accounts for the non-condensed bosons. Thus, 
/9o(x, x', i, 0) = LpQ{x.)ipo{x'){aQ{t)ao{0)) . In what follows we will omit the coordinate dependencies of po{x,x' ,t), 
and will call the correlator po{t) = (aJ(i)ao(O)) as the condensate (time) correlator. Accordingly, the condensate 
OFDM becomes po = p(){t = 0) . Apparently, po should be identified with the mean population Nq of the condensate. 

We assume that the equilibrium condensate can be well described by a macroscopic population of only one single- 
particle state ifQ. There are situations when this assumption is not valid. These are related to the low D geometries 
in which the condensate may exist as a quasi-condensate , or is thermally smeared over many states characterized 
by different winding of the phase jl^. In this paper we do not consider such situations. Thus, our analysis can be 
applied to the 3D case only. Furthermore, if there is the spin degeneracy, the so called fragmented condensate can 
be formed |2^. In this case, many states can also be macroscopically populated. Here we do not consider such a 
situation as well. 

Employing the standard definitions of the Heiscnberg operators as well as of the thermal mean, one finds 

po{t) = (4W«o(0)) = ^ p™e*[^-(^)-^-'(^-^)l*|(ao)™',™|', (2.2) 
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where the Boltzmann factor is Pm = Z ^{N, (3) exp{—PEm{N)), and Z{N, j3) denotes the canonical partition function 
as a function of the total number of bosons N and of the temperature T = 1/(3; the notation 

(ao)m',™ = {N - 1, m'|ao|m, iV), (2.3) 

has been introduced. Eqs. (|2.2| ), ( pT^ ) are exact. Below, we will represent po{t) as a cumulant expansion which is 
actually the 1/N expansion in the exponent, and we will find the leading term of this expansion. 

We note that several papers [2l|j24|l7|| have been devoted to developing an approach for treating iV bosons under 
the explicit condition of conservation of N . In this section, we introduce a natura l ph ase-space for the A^-conserving 



approaches. When N is conserved, it is convenient to express the matrix elements (2.3) in such a way that N becomes 



a parameter. First, we represent the exact eigenstate |m, N) of N bosons as an expansion 

\m,N)^ J2 CN,,N,..Xm,N)\No,Ni,N2,...), Nq ^ N - (Ni + N2 + ...), (2.4) 

Ni,N2,... 

in the Fock space \No, Ni, N2, ■■■) of the population numbers No, Ni, N2, ■■■ of some set of the single par ticle states 
tpoix.), (^i(x), (/32(x)..., respectively, given that </3o(x) is the only macroscopically populated state. Eq.( |2.4[ ) takes into 
account explicitly that the total number of bosons N is conserved; Cni,N2...{''^t-^) denotes the coefficients of the 
expansion. These coefficients form the Fock representation of the states with given N . We call Cni,N2...{''''^j ^) the 
projected states, and introduce a short notation |m, N) for them. Accordingly, the product of two projected states 
is defined as {N' ,m'\ra,N) = X^Wi at, N2 ('^^'i ^')^Afi,A'2---(™' Thus, for the same N, the product of 

any two projected states coincides with the product of the corresponding eigenstates ( ^.4| ), that is, {N,m'\m, N) = 

{N,m'\m, N). Therefore, the orthogonality condition {N,m'\m, N) = 6m'm holds. 

We note that, in general, no orthogonality exists between two projected states with different N. In other words, 
(iv7m|mV^') = J2ni.No....^*Ni,N2...(^' ^)^Ni,N2...{rn' , N') ^ for N' ^ N. This should be contrasted with the 



trivial orthogonality {N, m\m! , N') = for N' ^ N oi the eigenstates (2^) insured by the orthogonality of the two 

Fock subspaces with different N . In order to clarify the meaning of the projected states \m,N), it is straightforward to 
recall the radial part Rnxi''') of the full wavefunction Rn,L{T)YL^L^{9,(l)) of a particle moving in a central potential. 
Then, conserving TV is analogous to conserving angular momentum L. While Rn,L{f) and Rn',L{f) characterized 
by different radial numbers n are orthogonal for given L, no orthogonality exists between i?„ ^(r) and L'{f) for 

The eigen-problem H\m,N) — E„i{N)\m,N) is equivalent to finding the coefficients Cni,N2..X™i ■ Accordingly, 

this problem can equivalently be reformulated as H(N)\m, N) — Em(N)\m, N) in terms of the projected states, which 
are the eigenstates of the corresponding projected Hamiltonian H{N) with the same eigenergies Ejn{N) for given N. 
Resorting back to the analogy with the motion in a central potential, H{N) can be mnemonically viewed as the "radial" 
part of the total Hamiltonian. The projected states take care explicitly of the conservation of TV. As a matter of fact, the 
projected states form a basis for previously used iV-conserving approaches pl| , p^ , p^ . Employing them, it is possible, 

in principle, to construct explicitly the projected Hamiltonian H{N) — ^^771 \tt^t Em(N){N ,m\ in terms of the 
projected states. It should be noted that, despite being possible in principle, representing the projected states \m, N) in 
the original Fock basis by the coefficients Cni,N2,...{''^t impractical as long as the interparticle interaction is 

finite. The most productive way of constructing the projected states is in terms of the pair operators suggested in 
Rcfs. [ pljp^ . We, however, will first analyze general properties of the projected states without constructing them 
explicitly (Sec. III). For this purpose, any complete representation is suitable. Then, we will demonstrate the meaning 
of the projected states by employing an exactly solvable model of the bath of non- interacting oscillators (Sec.IVA). 
The role of interactions between the oscillators will be analyzed in sections IVB, IVC. We will also show (Sec.V) that 
the traditional hydrody nam ical approach |23,E4| provides a natural formalism for constructing the projected states. 



The matrix elements (2.3) can be expressed in terms of the overlaps of the projected states with N differing by 1. 



Indeed, a substitution of Eq.(2.4) into Eq.(|2.3D yields 



(ao)m',m — 

,N2...i'm,N), (2.5) 

Ni,N2.... 



where A'o is given in Eq.( |2.4| ). We consider a situation when the mean No of the condensate population A^o for given 
N is macroscopically large. As was shown by Giorgini et al. in Ref . , the relative mean square fluctuat ion SNq / Nq 
of A'o vanishes in the canonical ensemble for iV ^ 1. Thus, in the calculation of the matrix elements (pTsf) one can 
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perform an expansion with res pec t to SNq/Nq <^ 1 , and retain only the zeroth order term 
the accuracy o{6No/ Nq) Eq.(|2.5|) gives 



(ao)m',m — V ^0 Xr> 

where we have introduced the overlap of the projected states 

Xm'^m^ C^V,,Ar,...(TO',Af- l)CAr,,Ar,...(m,iV) = (iV^,m'|r7CiV) 

Ni,N2,... 



Nq, so that within 



(2.6) 



(2.7) 



with N differing by 1. In what follows we will assume that iVo ~ N. It is useful to mention some obvious properties 
of the overlap (2/7). Employing the orthonormality condition for the projected states with given N, one finds 



Xm' ,m"Xm,m" — <5ji 



(2.8) 



In Beliaev's work |2^ it has been shown that, for the exact ground state |0, N), one can write (ao)m',o — V No{Srn' ,o + 
o{l/N)). In other words, creation or destruction of one boson in the condensate does not lead to creation of excitations 
within the accuracy o {l/N ) (no external bath was considered in Ref. |^6|] ). Bel iaev 's result can be reformulated in 
terms of the overlap (|2.7| ) as xo.o = 1 + o{l/N). This immediately yields Eq.(2.2) as po{t) = iVo exp[i(i?o(iV) — 
Eo{N — l))t] + o(l) |26|] at T = 0, which implies no decay as long as N is fixed. 



At finite temperat ures T ^ 0, Beliaev's result is traditionally extended 
terms of the overlap (2.7), this is identical to 



as follows (ao)r 



Xr^ 



o(l/7V). 



,,. In 



(2.9) 



We introduce the notation 



— X]mP™l('^o)m,mP- Obviously, p'q^ ^ which plays an important role in the 
following analysis, is the diagonal part (with respect to the excitation label) of the condensate OPDM po = 



E™,™' Pml(ao) 



A^o- Employing Eq.( |2.6| ), it is convenient to represent p"^^ as 



Po — X , 



X = 



I / ^ 



(2.10) 



where we have introduced the mean overlap x- A direct consequence of (2.6) is 



X = l + o(l/iV). 



(2.11) 



The physical interpretation of Eq.( 2.11 ) is that, while entering or exiting the condensate, a boson does not significantly 

disturb the excitations reservoir. 

Employing Eqs.(2.6),( |2.9| ) in Eq.( |2.2| ), we obtain the condensate correlator ( ^.2[ ) as po{t) — p\^\t) +o(l), where 

we have introduced the diagonal part (with respect to the excitations) p'^\t) of the correlator ( |2.2| ). Given the 



condition (2_^), it is 



plf\t)^NoJ2Pme'''"'' +0(1), 



(2.12) 



^ dE,n(N) , , , 

/i„, = Era{N) - E^{N - 1) = ' + 0{l/N). 



Following the standard definition ||27||, the quantity /i„j introduced in Eq.(2.12) will be called the chemical potential 
of the mth exact A^-body eigenstate, so that the dephasing can be viewed as occurring due to the canonical ensemble 
fluctuations of the so defined chemical potential. 

Obviously, plf\t) « A^o ~ A^ for times short enough, so that no dephasing takes place. It is worth emphasizing that , 
while no decay is observed for each particular exponent representing the m-th eigenstate, the summation in Eq.(2.12) 
will normally produce a dephasing at times longer than some dephasing time tj_. We also mention that a statistical 
uncertainty inevitably present in the initial value of A^, namely: th e sho t noise (from realization to realization) , will 
result in the dephasing as well The evolution represented by ( 2.12| ) can be viewed as reversible dephasing in a 
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sense that no normal excitations are disturbed by such an evolution (see also in Ref. [^). We, however, note that 
this definition does not necessarily imply that the evolution can always be reversed in time prac tical ly. 

An alternative to the situation described above is that the mean overlap introduced in Eq.(2.10) becomes x — 
o{l/N) . It is important to emphasize that such a situation do es no t conflict with the existence of ODLRO. Let us 
discuss this. Apparently, if x = o[\/N), the diagonal part Po (2.1C), instead of being « N, becomes 



(d) 
Po 



No r = 0(1)- 



(2.13) 



On the other hand, the condensate OPDM po — Nq ^, Pm\Xm,m' = Nq ^ N due to the condition (p!q). In fact, 
the summation here runs essentially over m' ^ m, that is, the sum is collected from the amplitudes of the processes 
when one boson, which enters or exits the condensate, significantly disturbs the bath. Thus, it is con ceivable that the 
ODLRO, which can be expressed as po ~ N, can coexist with the situation represented by Eq.(2.13). We will refer to 



such a situation as the Orthogonality Catastrophe (OC) |50| in the system consisting of the Bose-Einstein condensate 
and the bath, which can be either extrinsic H or intrinsic 



We point out that, should the OC occur, the dynamics of the correlator (2.2) would be totally determined by Xm'.m 
with m ^ m', that is by the processes of creation and destruction as well as of scattering of the normal excitations. 
Consequently, the condensate correlator would exhibit a dissipative decay as long as the excitations have finite life-time 
due to the many-body interactions. 

The notion of the OC |^ was introduced by P.W. Anderson with respect to the Fermi-edge singularity where 
creation or annihilation of a single hole dramatically changes all states of the Fermi sea, so that the overlap of the 
states (projected on the Fermi-sea ) which differ in the number of holes by 1 is essentially zero in the thermodynamical 
limit. The OC in the bosonic system is defined above in terms of the mean square overlap x of the projected states 
which differ only in the total number of bosons TV by 1. Thus, in this case the role of the hole is played by 1 boson 
removed from (or added to) the condensate, and the role of the Fermi sea is played by the ensemble of the normal 
excitations (or by the external bath). 

Summarizing, if X ^ j the OC occurs, a nd t he condensate exhibits the irreversible decay. If x = 1 -t- o{l/N) , 
no OC occurs, and the condensate correlator (2.2) is essentially described by Eq.( ^.12| ), which, however, may exhibit 
a dephasing as a result of the thermal averaging. Below we will derive a general expression for x ior the case when 
N is macroscopically large. We will also obtain an expression for the correlator (2.2) in the main limit. 



III. OVERLAP AND THE SPECTRAL FUNCTION OF THE BATH 



If the total number of bosons N is conserved, the projected Hamiltonian H(N) can be constructed. In this 
Hamiltonian, N play s a rol e of t he parameter. Consequently, the mean overlap x between the projected states with 
different N (see Eqs.(2.7), (2.10)) can be found by employing the perturbation expansion with respect to 



H' = H{N) - H{N -I) 



dH{N) 
ON 



(3.1) 



Each additional derivative d.../dN in the exp ans ion (3.1) introduces a factor ~ 1/iV. Therefore, in the limit ^ 1, 
it is enough to co nsid er the first term only in (|3.l|) . Correspondingly, calculations of the overlap are based on the 1 /N 
expansion in Eq.(3.1) |^l|]. Below we follow the approach first introduced by Feynman and Vernon in Ref. for 
describing weak interaction with the bath. As pointed out in jH], it is sufficient to consider the effect of the bath in 
the second order with respect to this interaction, and the result must be considered as a first term of the expansion of 
the exponent. In our case the weakness is insured by the 1/N expansion. In other words, we look for the eigenstates 
|m, A^— 1) of H{N—1) assuming that the eigensta tes \m, N) and the eigenergies Em{N) oi H{N) are known. 
Accordingly, in the lowest order with respect to H' (^.l[) we find for (2.7) 



1 



Xm,m — 1 2 / ^ 



)mn I 



UJt, 



UJr, 



where the following notations 



{N, m\H'\n, N), uj^n = Er,{N) - Era{N), 



(3.2) 



(3.3) 



have been introduced, and the eigenstates \n,N) as well as the eigenenergies Em{N) are exact with respect to the 
Hamiltonian H{N). 
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As pointed out by Feynman and Vernon [Q, if there is an extremely large number of bath degrees of freedom, the 
exponentiation of ( |3.2| ) must be done after the averaging. This is a consequence of the central- limit theorem [0. 
Thus, averaging \xm,r. 



in Eq.(|3.2D over the ensemble and then exponentiating, we find the mean overlap 



X = exp 



where the spectral weight J{lo) is defined as 



J{uj) = y^^Pm\{5H),nn\'^5{uj 



(3.4) 



(3.5) 



and we have introduced the operator 5H whose matrix elements are 



{5H),nn — {H')mn — {H\ 



mm ^mn • 



(3.6) 



We note that the accuracy of the r esul t (3.4)- ( p.q ) is well controlled by the expansion (3.1). 

The OC occurs if the integral in (3.4) diverges at the low frequency limit w — > 0. The borderline for this divergence 
is the so called ohmic diss ipat ion |33y2[| characterized by the limiting behavior J(w) ~ (s = 1) as a; ^ 0. For 
s > 1, the integral in Eq.(3.4) is finite, which insures that no the OC occurs. In the case when the bath spectral 
function is ohmic (and subohmic) the OC will occur, so that the condensate correlator will irreversibly decay at large 
times. 

Let us assume that no OC occurs, that is, % « 1 in Eq.(3.4). Then, the correlator (2.2) is represented by Eq.(2.12). 
Its structure can be interpreted in terms of the thermodynamic fluctuations of the chemical potential Employing 
the central limit theorem, we find 



Noe 



o(l) 



where the mean value of /i„i and the dephasing time are given as 



(3.7) 



(3.8) 



and the definition of /i„i in Eq.(2.12), which can be written as — dEm{N) / dN + o{\ / N) , as well as the identity 
dE„i{N)/dN = {H')mm, with H' defined by ( ^T] ), have been employed. 

We note that the gaussian form (3/7) has its limitation time r', so that ( |3.7| ) is valid for t <C r'. It is importa nt to 
realize that r' ^ in the macroscopic system, so that wh en t is larger than by only few times the correlator (3.7) 
is essentially zero. Therefore, the gaussian approximation (3.7) is practically exact (if no spontaneous revivals of the 
kind Q are to be expected to occur at times ^ t^). We will discuss this r' in detail later. 

Note that the overlap {iA) can be interpreted in a different manner. Indeed, transforming Eqs. ( ^.5[ ), (3.4) into the 
time representation and introducing the spectral function G{t) in the time domain as the inverse Fourier transform 
of J{uj), one finds 



X = exp 



dt 



dt'G{t') , G{t) = {5H{t')5H{0)), 



(3.9) 



where the brackets denote the thermodynamic averaging over the canonical ensemble with given N; and 6H(t) — 
exp[iH{N)t]dH exp[—iH{N)t]. Thus, the OC is tightly connected with the long time behavior of the correlator 
G{t). If G{t) does not exhibit long-time correlations, the conjecture of the presence of the OC implies that G{t) 
contains the white noise part, that is G{t) — rd{t) + G'{t), where F > is some effective decay constant, and G'{t) 
stands for the part which satisfies dtG'{t) — 0. The decay constant F determines the time scale toe = 1/r 
on which the OC develops. Evidently, the validity of the traditional treatment of the bosonic system is limited by 
times t <ti toe (or by the excitation frequencies uj 3> toe) 113 ™der the assumption that the condensate correlator 
(4(t)ao(0)) =exp(-zp) ||. 

It is also possible to find the correlator ( ^.2[ ) as a 1/iV cumulant expansion, that is, 1/N expansion in the exponent. 
In the main 1/A'^ order, we find 



Poit)=Noe 



XW; X(i)=exp 



dti 



dt'Git') 



(3.10) 
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and, then, performed the exponentiation up to o{l/N); the 
.9|). Here we have introduced x{t) which can be interpreted as a 



where we have employed Eqs.(3.2), ( p.6[ ) in Eq. (|2. 
quantities JI, Td and G{t) are defined by Eqs.( |3.8[ ), 
time dependent overlap, so that x — x(i — oo). 

We emphasize that the gaussian factor exp(— t^/rj) in the form ( 3.1C ) is not due to the short time approximation of 
the correlator G{t). As the exp ressi ons (3.8) and ( ^^ ) indicate, the dephasing rate t'^^ is determined by the diagonal 
elements of the derivative H' (3A) of the projected Hamiltonian H{N) , while the correlator G{t) ( pT9| ) is built on 
essentially off-diagonal elements o f H' . 

The physical interpretation of ( 3.10| ) is the following. The condensate correlator po(i) is the ensemble mean of 
the correlators {N, rn\a\{t)aa{fi)\m^ N) in the given exact A^-body eigenstates. These partia l corr elators contain the 
non-decaying parts exp(j^mi)i characterized by the chemical potential /x^ (defined in Eq.(2.12)). The correlators 
{N,m\a\^{t)aQ((})\m^N) may also contain the decaying parts exp(— FmOi characterized by some decay constants 
Em > . In general, there is no direct relation between the parameters E^^ and /im: because is determined 
by the diagonal part of the operator H' (3.1), while F^ are given by the off-diagonal elements of this operator. 



The thermal averaging of the partial correlators should not change the exponential decay much, so that it can be 
char acterized by a single exponent ~ exp(— Ei) with some mean decay constant E, which can be related to G{t) in 
(3.9) as E « dtG{t) at times longer than some typical corre lation time of G{t) . It is important that the thermal 
averaging produces the gaussian time dependence in Eq.(3.10) due to t he fluctuations of in the non-decaying 
factors exp(i/imt). Thus, in general, one may expect the correlator (3.1C) to acquire the form 



-rt 



(3.11) 



at times longer than some typical correlation time of G{t) . We note that considering thermal (ensemble) fluctuations 
of /Ltm makes sense, if the dephasing time Td is shorter than (or at least comparable to) the OC time toe = 1/r. 

As we will see later, perturbative calculations yield toe N, while Td ~ \/A in the case when no e xternal bath 
is present. In this case, the OC becomes irrelevant for all practical purposes because po{t) of Eq.( 3.1C ) can be well 
approximated by the form (3.7) as long as t <^ toe, a-nd it will decay to almost zero at < i <C toe- Eurthermore, 
we will discuss that, in the isolated trap, the calculated toe turns out to be longer than an inverse of the typical 
value A of the matrix elements of the interaction part of the many-body Hamiltonian. This implies that the effect of 
the level repulsion should modify significantly the correlator G{t) at times t > . As will be shown, the effect of 
the level repulsion eliminates completely the OC in the isolated condensate, so that the evolution is dominated solely 
by the gaussian decay (3.7). 

Thus, while the necessary condition for irreversible decay of the condensate correlator ( ^.2[) is the OC represented 
as X = x(t = cx)) — !■ 0, a sufficient condition for the irreversible decay is toe ^ Td, provided the time scale on which 
the effect of the level repulsion may become significant is much longer than toe- This situation cannot be achieved 
in the isolated trap, and some external bath characterized by its own corresponding spectral function is required in 
order to create the OC. In the following discussion we will show that, as long as an infinite interacting bath contacts 
a finite condensate, the OC should always develop. 



IV. ORTHOGONALITY CATASTROPHE IN THE BOSONIC TRAP IN THE PRESENCE OF AN 

EXTERNAL BATH 

Eirst, we analyze the simplest model which exhibits the OC. In this model we neglect the normal excitations of 
the condensate, and take into account excitations of some external bath, represented by an infinite set of the non- 



interacting linear oscillators l32l 
constructed explicitly. We will a 



33| ]. In this case, the projected states as well as the spectral function J {to) can be 
so analyze the role of the interaction in the bath, and will demonstrate that the OC 



should occur in the interacting infinite bath even though no OC was found without such an interaction. 



A. Bath of non-interacting oscillators 

We neglect now excitations inside the trap and assume a presence of some external bath which, however, does not 
exchange particles with the trap. Eor illustrative purposes, the trap is represented by a single level containing all 
the bosons. The bath is described by an infinite set of linear non-interacting oscillators [ ^2[|3^ . This toy model is a 
simplified version of the model suggested by Anglin p5| . Thus, we have the Hamiltonian 

ffo = (eo + g)4«o + E[y + ^]' Q-Y.Ui^ + Y'^')^ (4.1) 

a a ^ ^ 
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where eg stands for the energy of a single level accommodating N bosons; the summation is performed over the 
non-interacting oscillators, with uja and ga, g'a being the frequency of the ath oscillator and the interaction constants, 
respectively. 



The advantage of the model (4.1) is that the eigenfunctions and the eigenenergies can be found explicitly. Fur- 
thermore, the pr ojected states can be constructed explicitly as well. The eigenstates and the eigenenergies of the 
Hamiltonian {iA) are 

|m, N) = |7V)|m>), |t7^>) = [] (g, - ga^-^N, il^), (4.2) 



7V2, (4.3) 



where \N) stands for the Fock state of N bosons occupying the level eo; = ^a{N) = ■y/cjJ'+g^lV, and we have 
explicitly indicated the dependence of on qa and fla- The projected states |m, iV), introduced in Sec. II, are 

actually the bath states which depend on as a parameter. The state |m, N) is represented here in terms of the 
oscillators eigenfunctions $„^(qo, fla), with m carrying the meaning of the set of all the integer quantum numbers 
{ria}- The projected Hamiltonian H{N) is, then, obtained from ( [4.l[ ) by simply replacing UqUq by N. Accordingly, 
the derivative ( |3.1[ ), which determines the spectral weight and the chemical potential becomes 

H' = H{N) -H{N-l) = e, + Y, (gaQa + ^ql) ■ (4.4) 

As discussed above in Sec. Ill, the diagonal elements of this derivative determine the chemical potential in the given 
eigenstate as = iH')mm + o{l/N) or 

= ^0 + E + ^) - il' (4-5) 

a " a 

where the label m now carries the meaning of the occupation numbers {ria} of the oscillators. The averaging of 
exp(— i^{„i t) ^ JJ^ ex.p{itg'^na/'2i}a) over {ua} produces the gaussian decay with the dephasing rate t^^ following 



from Eq.(3.8) as 



^-^"'^Eilt^a+^a+g]- (4.6) 

a ^ 

Here fia = (exp(/3a'a) — l)"-"^ stands for the mean thermal occupation of the a-th oscillator. It is important to note that 
in the canonical ensemble, the term ~ g'^ in Eq.( |4.5| ) does not produce any dephasing. If, however, TV fluctuates 
from realization to realization (shot noise) and does not change during the each realization, the dephasing rate 



(4.6) will acquire a contribution given by the variance AA^ of TV , so that the result ( |4.6| ) should be modified as 

The off-diagonal elements {H')mn define the spectral weight J(w) which may lead to irreversible decay. Which 
one of these two effects - either the reversible dephasing or the irreversible decay - dominates depends on the model 
parameters ga, g'a as well as on t he density of states of the bath 



Let us calculate the overlap (2.7) directly in order to see how the central limit arguments work for an infinite bath. 
This overlap is {N^l,m\r^) ^ JJa I dqaKA'la-9a^-^{N-l){N-l), 17„(TV-l))$„Jq,-5,l]-2(iV)iV, 1],(TV)) , 
and it can be found explicitly. There is, however, no need to calculate it exactly because the coefhcients ga, g'a must 
be scaled with the number of the effective degrees of freedom of the bath. Let us assume that the bath volume 
Vf, ^ oo is proportional to this number, and the volume occupied by N bosons is fixed. Then, the summation ... 
introduces a factor VJ,. In order to maintain the energy an extensive quantity, we impose g'^ and ga ^ 

This implies that the overlap, which can be represented as {N — 1, m|m, TV) — exp{^^ ln[/ dqa^n ila ~ 9a^a'^{N — 
l)(iV — 1), ^a{N — l))$„^(qa ~ 9a^a^{N)N, flaiN))]} , cau be expanded in the inverse powers of Vf, in the exponent. 
Accordingly, the only terms surviving the limit Vb — *■ oo are linear in g'^^ and quadratic in ga- This is a reiteration of 
the central-limit theorem |32|. Thus, the mean overlap ( ^.10| ) x acquires the form (3^), with the spectral weight 
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60 where (4.4) must be employed in the Hmit Vb oo. 
^ ~ l/VVb in the hmit Vh oo (if no shot 



Of course, the same expression ( |l7| ) follows from Eqs.(U), (| 

It is important that the dephasing rate (^) vanishes as 
noise is taken into account). The occurrence of the OC, which is expressed as the condition ( |2.13D , depends on the 
low frequency behavior of (4.7). If, e.g., J(0) ^ 0, the integral (3.4) will strongly diverge implying x ^ 0. Thus, 
the OC is equivalent to the irreversible dissipation for the considered model. Because, should the OC occur, its time 
toe ~ l/'^(0) does not contain any positive power of the bath volume VJ, , so that toe ^ Td in the limit considered. 
This corresponds to the irreversible loss of the phase memory on the time scale which is finite in the thermodynamical 
limit. A concrete realization of this situation, requiring an external bath, will be discussed elsewhere. 

If, however, the bath spectral function does not lead to the OC (e.g., J(w) ^ u)^ as w — > ), the time toe is 
formally infinite, and the only effect is the dephasing with tj. ~ a/14- We note that the last situation resembles most 
closely the situation in the isolated trap, where the role of the bath is played by the normal component confined in 
the same trap. 

Concluding this section, let us determine the validity of the gaussian approximation ( p.7| ) (or ( 3.10| )) wit h res pect t o 
the given model in the case Vb oo. Keeping in mind that ^ 1/Vt,, as discussed above, the form ( |3.7[ ) (or ( ^.10[ )) 
represents the first term in the exponent which is ~ l/Vt- Accordingly, it is valid as long as the next term is much 

It becomes of the order of one, when t > r' ~ Vt^^- Keeping in 
* '''dYb^^ ^ Td- Thus, as long as the system is macroscopic, the 



less than 1. The next term is Y.^ g[^t^ ^ i /Vb- 
mind that ^ \/V b (see (14.6D), we obtain that r' 



approximation (3.7) (or (|3.10|)) is practically exact, as it was discussed in Sec. Ill (see below Eq.(3.8)). 



B. The OC induced by an external infinite bath of interacting degrees of freedom. The strong OC. 
Now we consider the role of the interaction in the bath in producing the OC in the case when no OC existed without 



interaction. Referring to the previous model, let us assume that ( [4.7D is such that no OC occurs. What would happen, 
if the interaction between is turned on? Such a question is relevant for the confined condensate. As it will be 
shown below, the intrinsic bath of the non-interacting normal excitations in the trap does not produce the OC. Thus, 
it is important to understand if the OC would emerge when the interaction between the normal excitations is taken 
into account. 

Let us assume that the total Hamiltonian la H = Hq + Hint , where Hq is the "free" part, and Hint stands 



for some weak interaction. For example, the free part can be given by Hq (4.1), and the interaction is of the form 
Hint = J2abc9abcqaqbqc, with gabc being the corresponding interaction constants. In fact, a concrete form of the 
interaction does not matter much. We choose the third-power non-linearity as an example which resembles most 
closely the situation in the actual trap, where the most important interaction vertex contains three lines. 

The emergence of the O C in the presence of the interaction can be viewed as a n appearance of the exponential 



decay in the correlator Eq.(3.10) leading to the suppression of the mean overlap ( |3.9[ ). This happens if the interaction 
induces the decay of G{t) entering Eqs.( 3.ld| ),(|3.9|), so that, e.g., G{t) « exp(— 7t)G(0) with some 7 > 0. Then, 



the correlator po{t) (|3.1C| ) acquires the factor exp(— Fi) with F — G{0)/^ at times t :$> I/7. The gaussian factor 
exp(— t^/rj) in Eq.( p.lO| ), which is due to the ensemble fluctuations of the chemical potential, is not sensitive to 
presence of the weak interaction because the chemical potential may gain only small corrections in the c ase Q abc — > 0. 



Thus, the gaussian factor will not be significantly modified by the interaction. Accordingly, the form (3.11) should 
now describe the condensate correlator. 

We take another look at the emergence of the OC. In zeroth order with respect to Hint, there is usually a high 
degree of degeneracy between various multi-mode excitations. This degeneracy is represented by the condition cJmn = 
eI^\n) - E^^\n) ^ for some set n m, where the energy levels Ei°\N) are found for Hq. Here and below, 
the superscript denotes a quantity calculated in the zeroth order with respect to Hint- Let us first consider the 
result of the perturbation approach. If the interaction Hint is considered in the lowest order, the matrix elements 
iSH)ini (0 between some degenerate states may become finite with some statistically significant weight. Then, 
if the energy levels are taken unrenormalized, the summation close to the degeneracy point will make J(0) 7^ in 

r 
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Eq.(|3.5D. This yields the OC because the integral (3^) diverges as duj^^^ — > 00, implyng x — ^ 0. Accordingly, 



the time scale of the OC can be estimated as toe ~ J ^(0) = 7/^(0). 

It should be stressed, however, that the role of the interaction Hint is two-fold. On the one hand, it opens transitions 



in SH (3.6) between the degenerate states, and, on the other hand, it removes the degeneracy by introducing the level 
repulsion. While the first effect can be treated within the perturbation expansion, the second one is essentially non- 
perturbative. These two effects act in the "opposite directions". The first one tends to make J(0) ^ even though 
originally it might be that J(0) = 0. The level repulsion removes the degeneracy, and should result in J(0) = 0. Thus, 
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whether the OC does or does not occur depends on how the level repulsion modifies the spectral weight in the hmit 

The level repulsion is presently extensively discussed by the Random Matrix Theory (RMT) It has been 

shown that, due to the level repulsion, the probability P{uj) — {S{oj — oJmn)) to find two eigenenergies Em and 
En separated by the energy "distance" u! exhibits the universal behavior P{uj) ~ io^ for w — + 0, with s = 1,2,4 
depending on the general symmetry structure of the theory |^ . Recently, in Rcf. this feature has been explored 
to predict a dramatic change in the dissipation in a chaotic system in the presence of a magnetic field. 

We, however, note that it is possible to give simple physical arguments which justify the irrelevance of the level 
repulsion in the case when the external bath is infinite, while the condensate is finite. Indeed, the effect of the level 
repulsion becomes significant when the energy difference between two states is comparable (or less) with a typical 
value of the interaction matrix elements {Hint)mn linking these states. Then, some energy splitting A ~ |(^^mt)mn| 
will occur. Otherwise, no level repulsion should be taken into account. In an infinite bath characterized by the volume 
Vb — > oo, the largest matrix elements must be scaled as {Hint)mn Vb~" j where a > is some power determined from 
the requirement that the energy of the bath is an extensive quantity. Thus, A ~ VJ^"" —^ 0, and the level repulsion can 
be practically neglected because the physically relevant quantities, which describe the condensate, cannot depend on 
Vb oo . Specifically, the time toe, discussed in Sec. Ill, should be finite in this limit. In this case, no level repulsion 



should be taken into account, and the overlap (3.10) becomes strongly suppressed as x(i) ~ exp(—t/toc) ~^ 



as long as toe < t ^ A^'-. We will refer to such a situation as the strong OC. 

In the case of the strong OC, the spectral function can be modeled as J{lu) — Jq = const up to some high energy 
cut-off uJc- The constant Jq can be restored from the normalization condition J duj J{uj) ~ {{SH)'^) obtained from 
Eq.( |3.5D . Thus, Jq « a;~^((5i?)^), and, in accordance with the previous discussion, the time scale toe ~ l/>A)i during 
which the OC actually occurs can be found as 

^oc-j0^«A-\ (4.8) 

where the mean can be calculated in the zeroth order with respect to Hint ■ 

Thus, a contact of a finite condensate with an infinite interacting bath should practically always result in the strong 
OC. 

C. Absence of the strong OC in the condensate contacting no external bath. The Weak OC. 
Now let us consider the case of the intrinsic bath, that is, the normal component confined together with the 



condensate ( for a moment, we neglect the possibility that the dephasing effect may suppress the correlator (3.1C) 
much before the OC occurs). In this case, the level repulsion plays an important role 



In general, if the operator 5H entering Eq.( |3.5| ) and defined by Eqs.( |3.6| ), (3.1) is not simply ^ Hint, it is 
conceivable that SH links those degenerate states (with statistical significance) which are not mixed by Hint- In 
such case, these states remain degenerate (up to the corresponding order). Hence, the level repulsion would become 
insignificant, and the OC may emerge. We note, however, that such a case looks very artificial, when no external 



bath is present, because SH is simply the derivative (3.1) of the projected Hamiltonian. In the weakly interacting 



systems, the projected Hamiltonian can always be expanded in powers of N . This automatically implies that SH is 
statistically dependent on H{N). In other words, those states which are linked by SH with statistical significance 
are already mixed by the corresponding term in H(N), and their energies are thereby repelled from each other. 



As it will be seen in the next section, the operator SH, which determines the spectral weight (^), is actually 
given as N~^Hint, where Hint describes the interaction between the normal excitations. Thus, Hint is responsible for 
turning on of the transitions between the degenerate states as well as for the level repulsion. The above arguments, 
which justified the irrelevance of the level repulsion, when the external bath is present, are not valid anymore for the 
intrinsic bath because the typical splitting A introduced above (see Sec.IVB) is scaled now with N instead of the 

infinite volume of the external bath. The value of A ~ |(^^int)mTi| stems from the nature of interaction between the 
excitations. The processes responsible for this interaction are the pair scattering so that one member of the pair either 
enters or leaves the condensate (for references see as well as the following discussion in Sec.V). Thus, A ~ 
( or ~ 1/ \/]V for fixed density) . 



It is important to note that the OC time follows from Eq.(4S) as toe ~ because now SH = SHint/N is 
nothing else but the fluctuation of the interaction energy per one particle. Accordingly, the system is now in the limit 
toe 3> A~^, which is the opposite to the condition for the strong OC. This implies, that the matrix elements {SH)mn 



as well as the excitation energies uJmn in (3.5) experience strong renormalization caused by the level repulsion. 
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The above arguments imply that the mean overlap (3.4) (or ( |3.9| )) in the isolated trap may be sensitive only to the 
frequencies where the level repulsion takes place u> < A. This contrasts with the previous case of the external bath, in 
which the OC occurs as the strong OC, that is, far from the region whe re the level repulsion may become significant. 
Such a situation, when the OC occurs due to the integration in Eq.(3^) over the region lu < A, will be referred to as 
the weak OC. 

The actual time scale toe of such weak OC should always be too long on any practical experimental time s cale . 
For example, if the level repulsion modifies J{uj) from Jq to J{uj) ~ lu, the divergence of the integral in (3.4) 
is logarithmic. This immediately yields toe ~ exp(...A^) ^ A, implying that toe is practically infinity as long as 

^ 1. Nevertheless, let us discuss a possibility of the weak OC in the isolated trap. 

As a matter of fact, J{uj) in (3^) should approach zero not slower than ~ as w ^ 0, because the matrix elements 
{SH)mn in ( ^^ ) are strongly rcnormalizcd by the effect of the level repulsion. Indeed, if \n) is a subspace of the 

degenerate states corresponding to some energy Em\N) (with respect to the free Hamiltonian Hq), the new states 
and energies can be found by means of solving the corresponding secular problem which includes this subspace only. 
Thus, the matrix (-ffo)n,n' + {Hint)^^\i must be diagonalized in the subspace of \n) . The first part of this matrix 
is Em\N)Sn^ni duc to the degeneracy. Thus, the new spectrum and the new states follow from the diagonalization 
of iHint)^^\i alone. Accordingly, the renormalized {Hint)n,h' is diagonal at the degeneracy. In fact, close to the 
degeneracy, (-ffo)mn — Eln'Smn plays a role of the perturbation with respect to {Hint)mn. This implies that the 
renormalized {SH)mn ~ N'~^(Hi„t)mn ^ as ujmn 0. The matrix elements must be smooth in the parameters, 
therefore LUmli — > provide s a n atural scale for the off-diagonal values of the renormalized matrix elements {SH)^^- 
Hence, the spectral weight (3.5), which is ^ |(<5-ff)mnP gains an extra factor in the limit w — s- 0. This implies 
that the frequency integration in_(3.4) over the region uj < A does not yield any divergence. In Appendix B, we 
will estimate the mean overlap (3.45by employing the two-level approximation for treating the level repulsion. The 
conclusion is that this overlap X = 1 + o{1/^/N). 

It is worth noting that the situation with the irreversible decay of the condensate correlator is very special with 
respect to the picture of the decay of the normal excitations. Indeed, the decay time of the long living normal 
excitations scales with the typical size L of the system, which is i ~ N^^^ in 3D. It is much shorter than A^^ ^ Vn. 
Thus, the normal excitations decay long before the level repulsion may produce any effect. In contrast, the OC time 
toe calculated perturbatively with respect to Hint is toe -^i which is much longer than A~^ ^ ^/N. Thus, the 
effect of the level repulsion is crucial for the long time dynamics of the condensate. 

The above analysis shows that neither the strong OC nor the weak OC should be anticipated in the isolated trap. 
The results discussed above are not sensitive to particular details of the model. Nevertheless, below we will estimate 
explicitly the overlap x within the hydrodynamic approximation. 

To conclude this section, we emphasize that the question of the existence of the OC in a confined condensate is, as 
a matter of fact, irrelevant as long as no external bath is present. Indeed, the earliest time when the OC may become 
important is toe- This time turns out to be ~ (see above). On the other hand, the dephasing time ~ y/N (see 
Sec. VI). Therefore, the correlator ( |3.10 ) will decay long before the OC may become relevant. The mechanism of the 
decay of the condensate correlator should be distinguished from that for the normal excitations, which decay in 3D 
long before either the OC or the dephasing may take place. 



V. THE OVERLAP IN THE HYDRODYNAMIC APPROXIMATION 

The notion of the OC is well defined in the thermodynamic limit V oo and iV — > oo. Strictly speaking, the 
projected Hamiltonian H{N) should be constructed by employing the A^-conserving approaches |pl| , p^ , ITI . As was 
discussed above, the occurrence of the OC is determined by the low frequency (long time ) response of the excitation 
ensemble. Thus, in order to construct this response, it suffices to consider the Hydrodynamic Approximation (HA) 
which is known as correctly reproducing the behavior of the cloud at low frequencies and momenta p3| . The HA 
takes into account correctly the conservation of the total number of particles N as well as the vertex renormalization 
effects. A drawback of the HA is the occurrence of divergences at large frequencies. These are, however, non-physical 
and can be taken care of by introducing the upper cut-off at some energy lOc of the order of the chemical potential, and 
at some momentum kc of the order of the inverse healing length. On the contrary, the divergences at low momenta 



and frequencies are real and determine the peculiar physics of the condensate |23 2J]. Thus, in what follows we will 
employ the HA in order to construct the projected Hamiltonian H(N). 
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A. Dephasing due to non-interacting phonon gas 



We start with the classical Lagrangian 



L 



dt 2m 



t/(x) + /i 



(5.1) 



where the interaction term with g — 



> ( a is the scattering length) is taken in the contact form, with m being 



the atomic mass and U{x) standing for the trapping potential. The chemical potential /z is introduced to take into 
account the constraint 



In the HA, the ansatz 



* = v^e*' 



(5.2) 



(5.3) 



is employed in Eqs.(5.2), (5.1), where n, stand for the density and the phase, respectively. This results in the 
Hamiltonian 



H 



/dx{ 



(5.4) 



where the terms containing gradients of the density, which ar e no t relevant for the long wave dynamics, have been 
omitted. We note that, if viewed together with the constraint (5_^) in which N is just a c-number (not the operator), 
Eq.( ^.4[ ) represents the classical version of the projected Hamiltonian H{N). Keeping this in mind, in what follows 
we will use the notations H and H{N) interchangeably. 

In order to construct the non- interacting Hamiltonian Ha, the density is represented as 



n(x) = no(x) + n'{x), 



(5.5) 



where no(x) minimizes H (5.4), and saturates the constraint (5.2). In the Thomas-Fermi approximation (see in [p4|) 

noix) = {^I-U{x))/g, (5.6) 
with the boundary condition r7,o(x) = for U{x) > ii. Accordingly, the constraint (^.2|) , which becomes 



(ixno(x) = N, 



(5.7) 



yields the size of the condensate and the value of /i. Here we neglect the possibility of topological excitations in the 
condensa te. Thus, no winding of the phase needs to be taken into account. 

In Eg. (^7^), n'(x) denotes small perturbati ons (phonons) of the density which do not result in any change of N. 
Thus, J dxn'{x) — 0. Then, the Hamiltonian (5.4) takes the form H = Hq + Hint, where the quadratic part is 



Ho 



dx 



—noiVey + '-n 
2m 2 



provided ( p.6[ ) holds, and the interaction is defined as 

H,nt= [ dx^n'iVef 
J 2m 



(5.8) 



(5.9) 



Subsequent quantization of the Lagrangian (5.1) shows that the density n'(x) and the phase ^^(x) are the conjugate 
variables. The phonon contributions can be described in terms of the Fourier harmonics 0k and nj^ of 0(x) and n'(x), 
respectively ( 0(x) = exp(ikx)0k , n'(x) = exp(ikx)nj^), satisfying the commutation relation 



[6'k, n'^'] = -«(5k,k' 



(5.10) 
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Thus, one can express nj^. = i-^- The corresponding Hamihonian is actually the projected Hamiltonian H{N), in 
which TV plays the role of a parameter through ng given by Eqs.( |5.6| ), (5.7). The eigenstates of H{N) are the projected 

states \m,N), which can be expressed as some functionals $m({^k}, '^■0)- 

For sake of simplicity, we will consider a finite cubic box with the side L, and will set the trapping potential t o , 
with, e.g., periodic boundary conditions. Then, no(x) becomes just a constant uq — N/V, V — L^, and Hq in Eq.( [3.8| ) 
can be represented as 



Hn = 



k 



9 



92 



2 dOi^de^k 



h k rin 



2m 



2V ' 



(5.11) 



where 6]^ = 9'^-^. Taking into account Eqs. (|5.8| ), ( p.9[) , H' in (|3.1 

H' = dH{N)ldN 



IS given as 
gN 



k 



(5.12) 



Here the first term represents the global kinetic energy density (notice the factor V in the denominator) of the phonon 



gas. This term contributes to the dephasing rate ( |3.8| ) at finite temperatures. The second term, which is the 

mean interaction energy per one boson, does not contribute to the dephasing in the canonical ensemble. However, in 
case the evolution of the correlator (1.1) is considered for the case when the initial state is a mixtu re of states with 
different N (or the averaging is performed over different realizations of N) , the last term in Eq. ( 5.12| ) may contribute 
to the dephasing rate (3.8). We will discuss this later with respect to the effect of the shot noise inevitably present in 
the destructive measurements 

Let us neglect for a while the last terms in Eqs. ( 5.1l| ),( pT2 ). Then, comparison of Eq.( 5.12 ) with Eq.(14) indicates 
that the condensate and the normal component can be viewed as the case ga = 0, dj^ = 0, 5^ 7^ of the toy model 
considered in Sec.IVA. Specifically, in the absence of the interaction between phonons, the projected Hamiltonian 



( 5.11 ) obtained above in the HA is practically the projected Hamiltonian of the form ( 4J ) , with ga — ^ and = , 
= \/nog/m k. T he sp ec tral weight ( ^^ then follows as J{uj) w ^{g'a /^'i)n'i^{uJ - Via) ~ cj^T^/iV at large 
T, where Eqs.(^), ( 5.12D , ( p3.11 ) have been employed. Thus, no OC occurs due to the ideal phonon gas because 
J{ijj) ~ w^/TV, so that the integral in the overlap (3.4) does not diverge for w — > 0. Furthermore, due to J{uj) ~ l/A^, 
the overlap x = 1 + o{\/N). Consequently, the correlator ( p. 2D is given by the form (3.7), where the dephasing rate 
given by Eq.( [4.6| ) turns out to be tJ^ ~ T/ y/N (in Sec. VI, we will present accurate calculations of t^^)- The following 



discussion will be devoted to proving that the above simple picture is not altered practically by the interaction between 
phonons. 



B. Interacting phonon gas 



Now let us c onsid er the role of the interaction between th e ph onons. It is convenient to remove the dependence on 
no from Hq (5.11) and transfer it to the interaction part (5.9). We employ the scaling transformation 6'k = ^^^d'l^, 
with A — (no/rir)^'* where rir stands for so me ref erence density which will be kept constant with re spect to changing 
A^. As a result of this transformation, Hq (5.11) changes to X^Hq (we neglect the last term in Eq.(5.11)), where Hq 
has the form (5.11) in which 6'k is replaced by and uq is replaced by rir- This scaling transformation will change the 
interaction part (5.9) as Hint ^ X^^Hint- Thus, the total Hamiltonian can now be represented in the new variables 
as 



H 



TTl 



E 



1/2 



{H'g + i?-„(). 
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' 2 d9{^d9'_^ 



^ -3/4 .^2 
riQ \ iTi 



2m 



^k^^-k 



k,q 



9^...). 



(5.13) 
(5.14) 

(5.15) 



(k+q) 



In t he r epresentation ( p3.13[ )-( p).15 ), the matrix elements of the operator (5iJ, which determines the spectral weight in 
Eq.(|3.5D and which is given by Eqs. (3. 6), (p^), (|3.lD, become 
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{SH),nn — {iH[nt)mn - 5mn{H[n^)mm) ■ (5.16) 

It is important to notice the factor 1/A^ in this equation. 

The meaning of the scahng transformation introduced above must be discussed. The overlap (2.7) is calculated 



between two families of the eigenstates. Specifically, the change N ^ N — 1 produces the change of hq = N/V as 
no ^ no — 1/V. Thus, the overlap is to be found between the states ^m{{Oii},no — V~^) and <i>m({^k}, f^o)• We 
represent the overlap (pj) explicitly as an integral 



- J l?0k*:„({^k},no"^"')*m({^k},no) 



(5.17) 



over all the harmonics {^k}- The eigenstates ^m d^kj , np) can be expressed in terms of the eigenstates 
*m({^k}."o) of the new Hamiltonian + H[^^ (|5J^)-(|5J^) (where 9]^ = X'^O'^ and \ ^ [uo/nrf/^). 
Specifically, <i>m({^k}, "o) = A*^/^$^(A{0k}, no), where M stands for the dimension of D^k = Hk^^k • Then, 
Eq.( ^.17 ) takes the form 



c,^M/2 „ 

1 + y^^^kCia 



+ y){M,no-l^"^)<i>:„({M,^o), 



(5.18) 



where 5\ is the change of A as np cha nges by 1/V . In fact, it is enough to consider only the term lowest in 1/A^, so 
that 5\/\ = — (4iV)~^. Expanding ( 5.18 ) in 1/A^ and exponentiating the result |^ within the accuracy N~'^, we 
find the logarithm of the mean overlap as 



Inx = - 



M 



Bo{T) - Bi{T) - B2{T) 



(5.19) 



where the terms linear in 1/N cancel automatically, and the terms Bq{T), Bi{T), B2{T) represent the following 
contributions 



^ ■' m '' k,k' 



;({^k},no)0; 



92 



kCk' 



-KA{OvS.no\ 



(5.20) 



Bi(r) =4no 



4^) 



k 



C({^k},no)0k 



92 



dnodOk 



(5.21) 



B2iT)^~8nl 



1 



oe^ ^'*^{{ei,},no)-^^<^Um,no) 



(5.22) 



The first term in ( 5.1£ ) contains the formally divergent dimension M. This can be represented as M = J2k 
V J (fik 00. The divergence of the integral / d^k is a general high momenta divergence of the HA, and therefore 
it is non-phys ical. The integral must be cut off from above at the inverse healing length kc- Thus, M ^ N, and the 
first term in (5.19) behaves as ~ 1/A^ 0, and therefore it can be eliminated. 

The term Bo{T) represents a contribution which is finite in the absence of the interaction ( 5.15| ). Thus, it originates 
from the free pho nons as discussed above in Sec.VA. Accordingly, elementary gaussian calculations with the free 
Hamiltonian ( 5.14 ) yield 



BoiT) = i 



exp 



T ' 



k 



(exp( 



T ' 



(5.23) 



with hui\i ~ k being the spectrum of the normal excitations of the free Hamiltonian (5.14). The term ( 5.23 ) is ~ l/N, 
and it converges at small k in the 3D case. The interactio n ma y slightly change the "free" form ( 5.23 ) without 
changing the main scaling dependence on N. Thus, the term (5.20) can be eliminated from the consideration as well. 

The term (5.21) describes a contribution which is non-singular at the degeneracy points LUmn (see the discussion 
i n Sec . IV). The proof of this is presented in Appendix A. Thus, the term (5.21) contributes as ~ 1/iV in the overlap 
( ^.19 ). Consequently, we eliminate Bi{T) as well. 
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Finally, the term ( ^.22 ) is the one which is formally singular at the degeneracy poin ts. Thus, one can represent 
X = exp{— B2{T)) . It is possible to recognize (see the Appendix A) this form as Eq.(3^) with the spectral weight 



J(^) 



3 

4N 



J2 P" 



i.-^int)irin ' 



(5.24) 



is given by Eq.( 5.15| ). It is important to note the factor in Eq.(5.24). 

Above we have justified that the projected Hamiltonian H{N) as well as its space of states can be expressed in 



where i?,'„t 



such a way that the dependence on N is included in the interaction part H^^^ = H[^^{N) (5.15) only. T his im plies 
that, should the OC occur, solely the interaction would be responsible for this. It is worth noting that in (5.24), the 
matri x elem ents are taken between exact eigenstates of the "primed" Hamiltonian iJg + ( ^.13 ). 



Eq.(5.24) contains usual HA high momenta divergences. Such divergences, which are non-physical, can be eliminated 
by the following procedure. Let us replace x x/xo,0: where xo,o stands for the overlap of the projected ground 
states |0, iV — 1) and |0, N). It h as b een shown that such an overlap is essentially 1 (see the discussion below Eq.(2.8)). 
Thus, the renormalized overlap (3^) becomes 



X = exp 



duj 



(5.25) 



where the temperature dependence is shown explicitly in J{lu) — J{lu,T). Eqs.(5.25), (5.24) still have divergences at 
high momenta, which, however, vanish at T = 0. These divergences must be cut off from above at the momentum 
||,||. 

Let us first make naive estimates of toe from Eq. 



the form (4 
we obtain 



(4.S). That is, we ignore the level repulsion effect, and consider 
in the lowest order with respect to the perturbation theory, where 5H is taken from Eq.( ^.16" ). Thus, 



toe 



Pm\{Hl, 



(5.26) 



where Eq.(5.15) should be employed. The denominator of Eq.(5.26) is the mean square total interaction energy 
between the phonons. Thus, it is scaled as ~ as any extensive quantity, and Eq.(5.26) yields toe N. On the 
other hand, the typical matrix element A ~ {H-^^)mn (see discussion in Sec.IVC), where is given by Eq.( ^.15 ), 

is scaled as ~ 1/\/N. Thus, toe ^ A^^, and the strong OC cannot occur (see discussion in Sec.IVB). We note 
that the dephasing time Td ~ VN (see below Eq.( 5.12| )). Thus , reg ardless of whether the weak OC does or does not 
formally occur, the dephasing effect dominates the evolution ( 3.10|) (because Td <C toe in the limit N ^ 1). This 



implies that the question whether the OC actually occurs beco mes irrelevant because the correlator (3.1C) decays at 
times > Td which are much shorter than the time toe ( 5.26 ) of the possible OC. Nevertheless, in Appendix B, we 
have discussed the role of the level repulsion, and show that no the OC occurs. Furthermore, in Appendix B, we will 
evaluate the mean overlap x, and it will be shown that the only effect caused by the interaction between the phonons 
is that X = 1 + o{1/\/N), while in the non-interacting phonon gas x = 1 -I- o(l/7V) as it was concluded in Sec. VA. 

Thus, in the case of no extrinsic bath, the decoherence of the condensate correlator is essentially given by Eq.( [3.7[ ). 
Below we will estimate the dephasing time for the box. 



VI. THE DEPHASING TIME OF A CONFINED CONDENSATE 



Here we will calculate the dephasing rate of the condensate time-correlator as described by Eq.(3.7). The 
chemical potential /i™ — {H')mm , where H' is g iv en by E q.(b.l2 ), consists of two parts /i'-"-' and /i'^. The term 

and it is not associated with excitations. In the 
We note, however, that in the destructive 



(0) ^ 



gN/V results from the last term in Eq.(|5J), ( |5.1l| ), ( |5.1^ 



canonical ensemble, it does not produce any contribution to 
measurements N will necessarily fluctuate from the realization to the realization producing the shot noise 



(3., 



Therefore, the term fi^^' will contribute into the dephasing rate (3.8). To estimate this effect, we assume that 



the shot fluctuations are characterized by some known variance AA^. 
noise as 



(r'd 



(0)n-1 



gAN 
V2V' 



This yields the dephasing rate due to the shot 



(6.1) 
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It is important to emphasize that the shot noise, which is due to the uncertainty of A'', does not disrupt the phase 
coherence in the each reahzation. It should not be identified with the situation when N is not conserved due to the 
unwanted escape (or deposition) of particles from (to) the trap during the evolution. The latter effect is a source of 
t he w hite noise which erases any memory, and which leads to the exponential decay instead of the gaussian dephasing 

The term is due to the normal component. In order to obtain it, the spectrum of the elementary excitations 
ek{N) must be found. According to the preceding analysis, the processes of scattering (elastic and inelastic) of 
the exc itations do not contribute to the correlator (2_^) in the leading order iV ^ 1. Thus, in order to find 
in ( [3.8[ ), the interaction between the excitations can be essentially ignored. This will give the energy Em{N) as 
EmiN) = ek(-^)fik, where nk are the corresponding population factors of the elementary excitations. 

Strictly speaking, in order to find the effect (r^^ 7^ 0) for N finite, the iV-conserving approaches (2^j2^ should be 
employed. Above we have employed the HA, which is the long wave limit of these approaches. The HA has been 
sufficient for estimating the mean overlap in order to justif y th e absence of the OC. However, here we are interested 
in accurate calculations of the dephasing time Td from Eg. ([3.^). These calculations contain the integration over the 
momenta which are divergent on the high end, and which therefore cannot be accurately eliminated within the HA 
approach. Thus, the methods |^,^ is the only known alternative. However, we note that in the limit iV ^ 1, 
the differences vanish between the approaches pl| , p2| and the non-conserving N Bogolubov theory |Q . This implies 
that, if we employ the approaches pi| , p2| in evaluating the elementary excitation spectrum, we will gain corrections 
of the order 1/A^ compared with the Bogolubov method. Keeping in mind that the result in (3.8) is itself ~ l/N, the 
corrections to the spectrum will produce terms ^ which should be neglected. 

Thus, in what follows we will empl oy th e Bogolubov method and take ek as the standard Bogolubov spectrum 
(see [ ^8p4| ]). We note that the form (2.12) indicates that, while changing by 1, the number of excitations rik must 
be maintained unchanged. Thus, neglecting terms o{l/N'^) in Eq.( p.8[ ), we obtain 



(6.2) 



It is important to emphasize that the thermodynamic fluctuations, over which the averaging (...) is to be performed 
in Eq.(B.2), affect not only the term containing the population factors nt explicitly. Indeed, as has been discussed 
by Giorgini et. al |^^, the condensate fraction Nq exhibits anomalous fluctuations (SNq) . Consequently, these 
fluctuations must affect the spectrum of the normal excitatio ns £k depending on iVo. Below we will see that the 
contribution due to these fluctuations can be dominant in Eq.(6.2) for large N and not very low T. We represent 
fk = f'k{No) ~ e\^{No + SNq) as an expansion in a smallness quantity ^JJEW^ / N q ^ 1 around the mean value A^o of 
A^o- This gives 



deM) dekiNo) , dek{No)dSNo 



dN 



dN 



dNn dN 



dNdNo 



-o{5N^/N, 



(6.3) 



The fluctuation SNq is understood here as happening at fixed A^, thus one can employ dSNo/dN — SdNo/dN . 
We note that dNo/dN can be taken as dNo/dN — 1, because the corrections due to the interaction contribute to 
the higher order in the gas parame ter Both this expansion and n.k = rife + 5nk, where Uk stands for the mean of 
rik, should be substituted into Eq. (|6.2| ). In calculating the average, only the lowest terms in the fluctuations {{Sri^), 
{6Nq), etc.) must be considered. For large A^, the differences between canonical and grand canonical ensembles 
practically disapper with respect to the normal component. Thus, we employ ((rik — (rik))^) — Hk(rik + 1), where 



nk 



exp(f )-r 



(6.4) 



We neglect the terms of the order higher than 1/A^ in Eq.(6.2) (e.g., (SNon-k)). Close examination shows that each 
derivat ive of ek in Eq.(3.3) introduces an additional power of the gas parameter after the momentum summation 
in Eq.(6.2). Thus, we retain only the first two terms in the expansion ( |6.3| ). To proceed, we employ the approach 
psj in calculating the anomalous fluctuations. Thus, we choose the Bogolubov representation |3^ , |24| of the normal 
component as 



tp' = ^^(Ukttk - 
k 

2^k ■ 



(6.5) 
(6.6) 
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UkVk 



2ek' 



(6.7) 



where au (o;^) destroys (creates) a quasiparticle with the momentum k; ric ~ Nq/V stands for the condensate density, 
and the Bogolubov spectrum is 



21.2 



2m 



-, 1/2 



2 2 
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(6.8) 



Fluctuations of the condensate population SNq are related to the fluctuations of the normal component by means 
of the relation lEsI 



iVn 



Accordingly, employing Eqs.(3.5)-(6.7) and rik = a|j.Q:k in Eq.(6.2), we obtain 



where the notations are 



^1 -o2L^[—m-] "'^K + l)' 



2^V dNa 



X 2 



1 l^ dek{N^) _ 

'>\^ dNn 



.dNo 
' dN 



(6.9) 
(6.10) 

(6.11) 
(6.12) 



In the term ( 6.12 ) there is the factor {{SdNo/dN)^) which can be found from Eq.(x9). The change of by 1 does 
not affect the excitations accordingly to the above disc ussion. Hence, the quasiparticle operators in Eq.(^^) should 
be taken as independent of N. Thus, we find from (|6.9|) 



dN 



^ / Q Q Q \ 



(6.13) 



This representation should be employed in Eq.( |6.12) . As discussed in Ref. [g^, the main contribution into (SNq) 

comes from small momenta where u1 ~ vl ^ V Nq. Thus, Eq.( |6.13 ) yields dNo/dNo « No/2No. Accordingly, 

{{SdNo/dNof) ~ {SN^)/4nI Em ploying this as weU as the exphcit form of {SN§) for the box from Ref. fH], we 
obtain for the dephasing rate (6.1C) 



'1.8^1/2 6.0^2/3 

N ^ A2/3 



(6.14) 



where the limit T 3> /i is taken . The first term under the square root is due to the ensemble fluctuations of the 
population factors (see Eq.(6.11)). This contribution discussed by the authors in Ref. ||3^ is effectively gained at the 
momenta close to fee- These fluctuations do not practically contribute to the fluctuations of the condensate population 
A'o 1^^. The second term (see Eq.( 6.12 )) describes the contribution of the ensemble fluctuations of Nq determined 
by the low momenta collective modes [t25| . It can be seen, that the two terms under the square root can be dominant 
in different re gions of A^ and ^. 

The result (6.14) is obtained for the box. It is cle ar th at, apart from the numerical coefficients, a similar expression 
can be derived for the actual oscillator trap. Eq.( 6.14 ) can be employed for order of magnitude estimates for the 
conditions of the actual trap Q. At large A^, the first term behaves as ^ A^^^/^, and the second ~ N~'^/^, provided 
T/Tc = const and the condensate density no ~ A^2/5 Taking values typical for the experiment Aq = 5 • 10^, 

T « 50nK, a « 5 • 10" n 



10 cm we find ^ 



2 • 10-^ and Tc « S.S^i^no'VTO « SOOnAT. For these 
values, the second term under the square root in Eq. ( |6.14 ) is much smaller than the first one due to the temperature 
factor. In this case we find the dephasing time « Is. If, however, temperature is increased up to T/T^ = 0.5, 
the second term becomes dominant. In this case, the dephasing time becomes Td < 100ms, which is within the 
experimental range |9j. In this estimate we assumed no shot noise due to the uncertainty of the total value A^. If, 
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however, the variance AA^ in (3.1) is large enough, the total dephasing rate \J t^^"^ + (t^^"*)"^ can vary over a large 
range. 

We note, however, that the results obtained above should not be directly applied to the experiment 1^ in which 
the evolution of the relative phase of the two-component condensate has been studied. In a following publication we 
will modify the above analysis in order to apply it to the experiment |9t]. 



VII. DISCUSSION AND SUMMARY 



In this work we have applied the concept of the Orthogonality Catastrophe to a bosonic ensemble in order to treat 
temporal correlations of a confined Bose-Einstein condensate, consisting of a finite and fixed (albeit large) number of 
bosons. The occurrence of the OC turns out to be a prerequisite to the irreversible decay of such correlations. We 
stress that by saying "the temporal correlations of the condensate" we mean the correlator which is associated with 
the evolution of the global phase of a single condensate, or the relative phase between two condensates which do not 
exchange particles. The above analysis does not apply to the decay of the normal excitations. 

We have shown that the OC occurs if some infinite interacting bath contacts a finite condensate. On the contrary, 
if no external bath is present, the equilibrium evolution of the condensate correlator can be treated as though the 
normal component is not perturbed at all by such an evolution. The correlations are dominated by the reversible 
dephasing caused by the ensemble averaging over the realizations of the normal component. The effect of the level 
repulsion is shown to be significant in suppressing the irreversible decay at times much longer than the dephasing 
time. 

In general, relation between the three time scales - the time of reversible dephasing Td induced by the ensemble 
fluctuations of the chemical potential; the OC time toe (p^); and the time A^^ , during which the effect of 
the level repulsion takes place - determines the nature of the condensate decoherence. In the isolated trap, the case 
Td ^ toe and toe S> A~^ is realized. Thus, the decoherence is the reversible dephasing. In the presence of an 
external bath, toe <C A^^ , and the strong OC develops (the level repulsion is insignificant). In this case, the 
reversible dephasing can be observed if Td < toe ■ Ifj however, Td > toe : the decoherence of the equilibrium 
condensate correlator is completely irreversible. Realization of either case depends essentially on the nature of the 
external bath and its interaction with the condensate. The strong OC with short toe deserves special attention 
because in this case the nature of the low energy excitations may change dramatically p^ ]. 

Our results should be compared with those obtained by Graham in Refs. We note that in the phenomenological 
approach of Refs. |^,^, the classical field {^{x,t)) has been treated as a dynamical variable coupled to the reservoir 
of the normal excitations. A corresponding evolution equation with the white noise due to the normal component 
has been postulated. The conclusion has been given that a major reason for the decay of the field {'i'{x,t)) is 
the irreversibility developing at times longer than some decay time containing a factor larger than y/N (see, e.g., 
Eqs.(7.48)- (7.50) of Ref. ||]). In this regard we note that, the assumption that {^{x,t)) ^ is a classical variable 
which satisfies classical equations of motion inevitably violates the conservation of N on the level of an uncertainty 
~ ^/N. Thus, the classical Bogolubov treatment is accurate within l/\fN . This does not introduce any problem unless 
the long time limit is considered (times ~ and longer) fl|l5|-|l7|. Thus, we consider the results of Refs. as 
based on an unjustified use of the Gross-Pitaevskii equation as well as of the Bogolubov approximation in the domain 
where these approaches should not be employed. In contrast, our treatment takes into account the conservation of N 
explicitly, and the Bogolubov approximation (in Sec. VI) was employed only for calculating the dephasing rate which 
itself has a prefactor \/\/N, so that the corrections ^ 1/N due to the Bogolubov approximation are insignificant. 

In future work we will extend our results to the case of the specific situation realized by the JILA group |^ . We will 
also discuss how the reversibility of the dephasing can be revealed in the echo-type experiments |3^]. It is also worth 
considering the possibility of creating a situation in which a confined equilibrium condensate contacts some external 
bath inducing fast irreversible decay of the condensate correlator . 
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APPENDIX A: TERMS Bi{T) AND B2(T) IN EQS. (5. 21), (5.22) 



The terms Bi{T) (|5.2l| ), and B2{T) ( p.22| ) can be written as 

^ ' m k 

and 



(Al) 



(A2) 



because f&^df^k}, ?^o) ^^re the projected states \N,m). In Eq.(|Al[) wc employed the standard representation 



d6[. 



(flk - ak)\/^ 



Zk 



mg 



(A3) 



in terms of the phonon annihilation and creation operators Ok and a|., respectively. These operators diagonalize the 
free part (5.14) of the projected Hamiltonian ( ^.13| )-( 5.15| ), so that Hq — X)k with the bare spectrum defined 



as Efe — Cok, Co = yh nrg/m. 

The matrix elements {N,m'\-i^\N,m) can be found by differentiating the relation H{N)\N,m) — Ern{N)\N,m) 
with respect to uq (note that no — N/V is the only combination which contains the dependence on N). Then, we 
find 



— 9 — 
{N, m'\- — \N, m) 
ono 



Ano E„,{N) ~ E„AN) 



(A4) 



for m ^ vn! and zero otherwise. In here we employed an explicit form ( 5.15 ). We note that (A4) is exact. Then, (Al) 
can be represented as 



i?i(T) = 



2(4A^)^ 



E^{N)-E^,{Ny 



(A5) 



The expression ( |A5| ) can be evaluated by employing the perturbation theory with respect to H[^^ (5.15). We note 

that in the lowest order, Bi{T) = beca use the zeroth order intermediate states |iV, m') differ from the state |iV, m) 
by two phonons, while {H'^^^^)„^'m ( 5.15 ) links those states which differ by either one or three phonons. The higher 
order terms do make Bi{T) ^ 0. However, the dominator Em{N) — Em'{N) in t his expression may lead only to the 
first order pole which does not lead to the divergence ~ l/ui'^ in the overlap (3.4). Thus, the term i?i(T) ^ 1/N. 

The term B2(T) in (|A2|) can be rewritten by employing the relation {N,m\£^\N,m) = -(g^(iV, m|)(g|^|A^, m)) 

as well as Eq.(A4). Finally, we find 



B2{T) 



3 

47V 



E P" 



intJm'm] 



■{E,n{N)-E„AN))^- 
This expression coincides with i J dio J {u) / oj'^ where J(w) is given by Eq.( pS.24 ). 



(A6) 



APPENDIX B: LEVEL REPULSION AND ABSENCE OF THE OC IN AN ISOLATED TRAP 



i„t;mn|, where the index ^^'^ indicates that the matrix elements 
At high T one can estimate 0^ ^ \/ rriT/h^k'^nr and 



We evaluate the typical matrix element A w | (i/,-, 
are taken in zeroth order with respect to i?j'„j ( ^.15 ). At high T one can estimate 6^ ^ \JmT/h?'k'^ 
d.../d9k ^ \fTjg from (5.14). Accordingly, we arrive at the estimate of A from Eqs. ( ^.13 )-( 5.15| ) as 

1/2 



T 
N 



(Bl) 
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This value sets up an energy scale below which the level repulsion becomes significant. Here the factor indicates 
that no the strong OC occurs. Indeed, as discussed in Sec.IVB, the necessary condition for the strong OC is that 
the OC time toe ( 5.26| ) (calculated by means of the perturbation expansion) is shorter than A^^. Eg. ( 5.26 ) yields 
toe ^ N, implying that for large enough N, the opposite limit holds toe ^ A"-'^. Hence, the OC in the isolated 
trap may occur only as the weak OC (see Sec.IVC). 

The relevant frequency scale, at which the effect of the level repulsion dominates, is given by 



LO < A. 



(B2) 



At such frequencies the level repulsion effect strongly modifies the spectral weight (5.24) with respect to the pertur- 
bation expansion. As discussed in Sec.IVC, a strong renormalization of the matrix elements takes place. This effect 

is non-perturbative. This can be understood within the effective two-states model, in which a pair of the degenerate 

- (0) - (0) P — . 

states \m,N) and \n,N) defined with respect to the free Hamiltonian Hq ( p. 141 ) is considered independently 

from the rest of the states. 

In order to justify this two-state approximation we note that the degeneracy contributing to the spectral weight 
J{lo) (5.24) at w — > in the lowest order of the perturbation expansion is associated with the processes of the decay 
and the recombination of two phonons. Given an initial phonon with some momentum k, the final states contain two 
phonons with the momenta p and q, which obey the conservation condition k = p -I- q. The subspace of the final 
states, which may become degenerate with the initial state, is characterized by different p and q. Only those final 
states, which are close enough in energy to the initial state, participate practically in the renormalization cause d b y 
the level repulsion. A natural criterion for the closeness is given by the typical value of the matrix element A ( |Bl| ). 
Let us assume that some pair p and q satisfies this condition. Then, in a finite system, another closest final state 
will be characterized by the momenta differed by the amount ^ l/L 0, where L is a typical size of the system. 
Accordingly, their energies should also differ by some amount ~ l/L 0. The value of A (Bl) is scaled, however, 
as ~ l/L^/'^ <^ l/L. Thus, for large enough system, such a closest state turns out to be out of the range of the level 
repulsion. These considerations justify that only pairs of the degenerate states participate practically in the level 
repulsion. 

.(0) „ . ... —.,(0) , , —,(0) 



Then, the matrix element {H[^^)., 



which is defined with respect to the bare states |m, TV) ' and |n,A^) 
< A). The new states |m, TV) and \n,N) are the linear combinations of the 



mixes these states (as long as 

- — - (0) - — - (0) 

bare states \m,N) and \n,N) . Then, the standard treatment of the two-level system close to the degeneracy 
shows that the renormalized matrix element {H[^^)mn and the energy difference LUmn become 



I,/") I 



iH' 

\ int/mn ' 



± 



(B3) 



for |wmi| <C A. This feature alone leads to a strong suppression of J{uj) (5.24) in the 



Thus, {m^,) 
range (p^). 

Let us estimate the mean overlap x ( p3.25 ) where the Jiuj) is given by Eq.(5.24). We introduce a distribution 

function of the bare matrix elements | (Hmt)mn | and of the excitation energies \ijOmn \ '■ 



(B4) 



so that the spectral weight (5.24) can be represented as 



J{uj) = 



AN 



dM 

^0 



duJoP{M,ujo)'^M^S{uj - JnP+ujI), 



(B5) 



where we have employed Eq. (|B3| ). 

We assume that, if no level repulsion is taken into ac count (that is, the matrix elements as well as the excitation 
energies are taken unrenormalized), the spectral weight ( 5.24 ) can be well taken as a con stant. This would imply the 
OC. A simplest form of P(Af, wq), which will produce such a spectral weight in Eq.(Bf:), if the level repulsion were 
ignored, is 



P(M,wo) = ^'oexp(- 



(B6) 
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where Pq is some normalization constant which can be found from the sum rule 



(B7) 



following from Eq.(5.24). The form (B6) indicates that the distribution of the matrix elements and the excitation 
frequencies is uniform, as long as the matrix elements are less than some typical element A (Bl) and the excitation 
frequencies are below the high energy cut-off lUc- Employing (B6), the unrenormalized spectral weight acquires the 
form 



dM 



dujQ e 



M'S{cj -u;o) = 2 [ — ] A'Poe 



(B8) 



It satisfies J(0) ^ 0. In order to obtain Po, it is enough to employ this unrenormalized form in Eq.(p7p, because 



the renormalization ( |B5| ) bec omes important only for w < A , where A in Eq.(Bl) is small as ~ \/\N. Then, a 
substitution of (B8) into (B7) yields 



1 



2tj,A3 



(B9) 



This can now be employed in Eq.(B6), and then in Eq.(B5). This yields the renormalized J{lo) (BE) as 



2wc 



(BIO) 



where the condition <C cj^ is taken into account. Thus, J{oj) ^ oj^ /N for cj <C A , and J{uj) « const ^ 1/N for 
A <C w ^ Uc- Such a beh avior i mplie s that the OC is eliminated due to the level repulsion. This can be seen directly 
after employing Eq.(BlO) in Eq.(5.25) for the mean overlap. We find 



Inx 



28 N^uJcA 



(Bll) 



A quick glance on Eg. ( Bll ) shows that Inx ^ — l/viV 0. Indeed, {{H[^^) ) is the mean square interaction energy 



~ N. The value of A is ~ l/vN (Bl). Thus, due to the iV^ in the denominator, the total factor follows. 
Accordingly, x = 1 + o{1/^/N) implying no the OC. 
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